Star product and ordered moments of photon 
creation and annihilation operators 



S N Filippov 1 and V I Man'ko 12 

1 Moscow Institute of Physics and Technology, Dolgoprudnyi, Moscow Region, Russia 

2 P. N. Lebedev Physical Institute, Russian Academy of Sciences, Moscow, Russia 

E-mail: sergey . f ilippovOphystech . edu, mankoOsci . lebedev . ru 

Abstract. We develop a star-product scheme of symbols defined by the normally 
ordered powers of the creation and annihilation photon operators, (a') m a n . The 
corresponding phase space is a two-dimensional lattice with nodes (m, n) given by 
pairs of nonnegative integers. The star-product kernel of symbols on the lattice 
and intertwining kernels to other schemes are found in explicit form. Analysis of 
peculiar properties of the star-product kernel results in new sum relations for factorials. 
Advantages of the developed star-product scheme for describing dynamics of quantum 
systems are discussed and time evolution equations in terms of the ordered moments 
are derived. 
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1. Introduction 

This article concerns quantum states of one-mode radiation field and their time 
development by using normally ordered moments ((a^) m a n ) of photon creation and 
annihilation operators, a) and a. Ordering of creation and annihilation operators was 
studied extensively in the 1960s in connection with the quasiprobability distributions 
(see, e.g., [IH3]). Thus, the Husimi-Kano function Q(a) [HE] is closely related 
with the normally ordered density operator p = J2 m n 

namely, Q(a) = 

^ mn Pmn(ct*) m a n , where a* denotes the complex conjugate of a. Similarly, the 
Wigner ty-quasidistribution [6] and the Sudarshan-Glauber P-quasidistribution [TUB] 
are expressed through symmetrically and antinormally ordered representations of the 
density operator, respectively (see, e.g., [SI [10]). The quasidistributions on phase 
space do not limit to Q, W, and P-functions and are reviewed in several papers (see, 
e.g., [101412] ) . The recent detailed review of the phase-space approach is presented by 
Vourdas [13]. An advantage of the phase-space description of quantum states is that one 
deals with functions (C i— > K) instead of operators. As equations for quasidistributions 
do not involve any operator, they are sometimes easier to handle than the Schrodinger 
or von Neumann equations [10]. Moreover, in the phase-space formalism, quantum 
phenomena are known to be interpreted in the classical-like manner [H]. 

Being a rather good alternative to the density operator, the quasidistribution 
functions Q, W, and P are all equivalent in the sense that they all contain the thorough 
physically meaningful information about the state pH]. As the quasidistributions are 
equivalent, the question arises itself: why do not we use only one of them for all the 
problems? The answer lies in actual applications of the quasidistributions. For example, 
using the classical interpretation of the Wigner jy-function in collision problems, it is 
possible to make some reasonable approximations and obtain an appropriate but still 
accurate solution to the problem without requiring excessive computer time and expense. 
On the other hand, the nonnegativity and smoothness of the Husimi-Kano Q-function 
make it advantageous for the analysis of classically chaotic nonlinear systems pi)]. Also, 
the Q-function of the radiation field can be measured at optical frequencies via an eight- 
port homodyne detection scheme (see, e.g., [9]). Finally, the negativity of the Sudarshan- 
Glauber P-function is a commonly accepted criterion of the state's non-classicality, 
however, to observe the fact P(a) < is not that easy (see, e.g., [IH1[T6]). We can 
draw a conclusion that although all the quasidistribution functions are equivalent, the 
different functions exhibit different properties and obey different dynamical equations. 
Those differences specify the most advantageous representation for a particular problem. 

The set of normally ordered moments ((a^) m a n ) can also be treated as a 
quasiprobability function of two nonnegative integers, i.e. / : Z + x Z + h-> C. An 
interest to the ordered moments rose in the 1990s and was encouraged by the advances 
of measuring quantum states of light. Similar to the conventional quasidistributions, 
the moments {(a^) m a n ) thoroughly determine the density operator p and the explicit 
relation is found in the papers [TT1 - [2~T] . Despite the moments ((a^) m a n ) are but one of 
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plenty equivalent prescriptions to deal with quantum states (see, e.g., [13]), there is a 
particular problem where the moments play an indispensable role: a measurement of 
itinerant microwave quantum states. 

The matter is that at optical wavelengths there is no need to calculate the normally 
ordered moments and then use them to find the density operator. The output of 
the homodyne detection scheme is rotated quadrature distributions w(X, 9), X e R, 
6 G [0, 7r], also referred to as optical tomogram (see, e.g., [22H21] and the review [25]). 
There exists an explicit formula for reconstructing the density operator p in terms of 
the optical tomogram (see, e.g., [26T428] ). The normally ordered moments are easily 
expressed through the measured optical tomogram as well [291130] . 

At microwave frequencies, conversely, to measure the rotated quadrature 
distributions w(X, 6) is a challenge [31] because of the strong thermal noise added 
and the absence of single-photon detectors in such spectral region. However, it has 
been reported recently how to extract the lower-order moments {(a^) m a n ) of microwave 
quantum states themselves and exclude the noise contribution to experimental data [321 - 
l3l] . For details we refer the reader to the paper [35J. Thus, in the microwave domain, 
the lower-order moments are experimentally determined and contain the primary 
information about a quantum state. 

Given only normally ordered moments of the microwave radiation field, it is 
reasonable to associate a quantum state with the measurable quasiprobability function 
f(m,n) = ((a)) m a n ) of two nonnegative integers and consider the time evolution 
f(m,n;t) governed by a Hamiltonian H in presence of decoherence processes. This 
motivates us to follow a star-product approach [361412] and develop the particular 
star-product scheme on a two-dimensional lattice (m,n), where any operator A is 
associated with a symbol /_a(^ ; n ) — Tr[(a^) m a n A] and the product AB of two operators 
corresponds to a star product (also known as twisted product) of the corresponding 
symbols, i.e. /ab(wi, n) = n). Hence, one can associate the density operator 

p of a quantum state with the normally ordered moments ((a)) m a n ) = f p {m, n) and the 
Hamiltonian H with its symbol Tx[(a)) m a n H\ = n). Then the unitary evolution 

of a quantum state can be easily written in terms of the introduced symbols as 
df (m Ti) % 

— JL Qf — = —fi[fH*f P -fp*fH](m,n), m,n = 0,l,..., (1) 

where h~ is the Planck constant. 

Formula ([1]) is nothing else but the evolution equation for the normally ordered 
moments. Therefore, the time development of a quantum state is expressed through the 
measurable characteristics, which resembles the expectation value approach [131111] and 
the tomographic-probability approach [28lll5j|l6] t° quantum mechanics. According to 
the latter one, the measurable tomographic probability is a primary object to describe 
quantum states. The evolution equations for tomograms were found, e.g., in [UH19]. 

A kernel that determines the star product * in formula (CQ) is not known in the 
literature and is expressed through factorials in Section 131 Since the kernel of any star 
product scheme is to satisfy a nontrivial sum relation, we derive a new relation on 
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factorials in passing. 

Let us note that a unitary evolution of any quasidistribution function can be written 
in the form of Eq. (CQ) (see, e.g., \±0\) but the specificity of our particular case is that 
the symbol /#(m, n) is not determined even for the harmonic oscillator Hamiltonian 
H = ftio{a)a + 1/2) because the trace is diverging. Nevertheless, in Section H] we show 
that the symbols [/# * f p ] (m, n) and [f p * f H ] (m, n) can still be determined and Eq. ([I]) 
takes the form of a difference equation for f p (m, n). 

The aim of this paper is to develop the star-product scheme of symbols given by 
the normally ordered moments, to explore properties of the star-product kernel, and 
to derive unknown evolution equations for the moments via the star-product approach. 
It is worth pointing out that we study the star-product scheme not of functions on 
the conventional phase space (q,p) but of functions on the two-dimensional lattice 
(m, n), where m,n are nonnegative integers that determine the moments {(a^) m a n ) . 
This peculiarity provides a nontrivial change of variables and, as we will see, yields the 
difference quantum basic equations in contrast to the partial differential equations for 
the quasidistributions Q(a), W(a), and P(a). 

The article is organized as follows. 

In Section [2j we show how to treat operators as symbols defined by the 
normally ordered powers of creation and annihilation operators. Also, we revisit the 
reconstruction of operators given by their symbols, i.e. the problem of moments. In 
Section [3j we develop the star-product formalism (JlJEE] °f the introduced symbols. 
In particular, we find the star-product kernel and intertwining kernels. We calculate 
those kernels in the explicit form and use their particular properties to derive new 
sum relations involving factorials. In Section HI we apply the developed star-product 
scheme to derive the time evolution of moments. Unitary and non-unitary evolutions 
are considered. Conclusions are presented in Section |5j 

2. Quantization scheme: symbols of operators 

Unless otherwise stated, we assume that the expression Tr[-] is well defined. Thus, we 
intensionally avoid discussions of the convergence problems and focus our attention on 
the recipe to deal with symbols of operators. 

Consider an operator A acting on the same Hilbert space as the density operator 
p. Assuming the existence of the trace 

f A {m,n) = Tr[(tf) m a n A] = Tr[U\m,n)A], (2) 

we will refer to the function n) of two discrete variables m, n — 0, 1, . . . as symbol 

of the operator A. For the sake of convenience, we have introduced the dequantizer 
operator U(m,n) = (a)) n a m . Symbols /^(m, n) are defined on the two-dimensional 
lattice (m, n) (Fig. An example of the symbol f p (m, n) is illustrated for a squeezed 
vacuum state p in Fig. [lb. 
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Figure 1. (a) Phase space is a lattice whose nodes are labelled by two nonnegative 
integers (m,n). Symbol / p (m, n) is a complex valued function on the lattice, (b) 
Vectors (Re/ P (m, n)/m!n!; Im/ p (m, n)/m\n\) represent complex values of the symbol 
/ p (m, n) of the squeezed vacuum pure state p with squeezing parameter r = 1.3e l5Tr / 4 . 



The set {/^(m, n)}^ n=0 is known to be informationally complete [PTH2T] . In other 
words, the operator A can be reconstructed as follows: 



A— fA(m,n)D(m,n) 



m,n=0 



v ' mini ^ (m + j)l(n + j)l v ' 

j=-{m,n} 



(3) 
(4) 



where {m,n} = min(m,n) and the operator D(m,n) is called quantizer. 

Using the Fock state representation of powers of creation and annihilation operators 
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the dequantizer and quantizer can be rewritten in the form 

^„ ) = f#ipl|„ +i)(m+k | 



D(m, n) 
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k=0 
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The dequantizer and quantizer are shown to be orthogonal in the sense of trace 
operation [T7H2T] 

Ti[U\m,n)D(m',n')] = 6 m , m 'S n ,n' (8) 
where Sij is the conventional Kronecker delta-symbol. 
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Let us now check that Eqs. (jnjHCH) define an informationally complete scheme. It 
is shown in the papers [50H52] that a quantization scheme is informationally complete 
(tomographic) if and only if 

oo 

\\D(m,n)))((U(m,n)\\=i, (9) 

m,n=0 

where X is the identity super-operator, \\D(m,n))) and \\U(m,n))) are vectors 
constructed from the quantizer and the dequantizer by the procedure of representing 
matrices as vectors. In our case we have 



\\U(m,n))) = If J \n + k)®\m + k), (10) 



fc=0 
{m,n} 



||D(m,n))) = 2^ — T7r = ^T7 = TT- 
Direct calculation yields 



OO 



(p|<g>(g| ll- D ( TO ' n )))(( c/ ( m ' n )H i r )®l s ) 



vm,n=0 




r\s\ -A f-l) n -f 



(r — n)!(n — »)! 

n=p v 7 v ' 



r!s! (1 - 1)*-" _ r!s! 

q\p\ [r — p)\ y q\p\ 

that is the relation (jUJ) holds true indeed. 
3. Star product 

By definition, a star product of symbols J'a and fs is nothing else but the symbol 
fAB{ m , n ) of the product of two operators A and B, namely, 

(/a* f B )(m,n) = fAB(m,n) 

oo 

= fA(m\n')fB(m'',n'')K(m,n)m\n';m'',n''), (13) 

m',n',m",n"=0 

where the kernel K is expressed in terms of the dequantizer and quantizer operators as 
follows: 

K(m, n; m', n; m" , n") = Tr[U\m, n)D(m', n')D(m", n")\. (14) 

From definition f[T3"j) it follows immediately that the star product is associative and the 
star-product kernel satisfies a nontrivial relation 

oo 

K(m, n- k, I; rri", ri")K(k, I; rri, n'\ m", n") 

o 

oo 

Y K{m,n;m',n']k,l)K(k,l]m",n";m"',n"'), (15) 



k,l=0 

oo 



k,l=0 
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Figure 2. The kernel K(m,n;m' ,n';m" ,n") is nonzero if the points (m,n), (m',n'), 
and (m",n") determine a parallelogram such that its forth vertex, that is opposite to 
(to, n), belongs to a bisecting line. Also the conditions m < m", n < n' are to be met. 



which is a consequence of the relation f A * f B * f c = (f A * f B ) * f c = f A * (f B * f c ). 

Substituting (JBJ and (J7j) for the dequantizer and quantizer in (JHJ, after some 
algebra, we obtain the star-product kernel in the explicit form 

K{m } n;m ,n ;m ,n ) 

{-l) m '- n "{m' + m" -m)\ 

71 77u 7/ \T7 ~i rrUm+n'+n" ,n+m'+m" 

m'ln "\{m" — m)\[n' — n)\ 
_ (-l) m -""r(m / + m" -m + l)5 m+nl+n „ , ra+m / +m » 

" ^(m / + l)^(n / ' + l)^(m"-m + l)^(r^ , - n + 1) ' ^ ^ 

where V is the conventional Euler gamma function. 

Using the property of the Kronecker delta-symbol, the kernel can also be rewritten 
as follows: 

K(m, n; m' , n'\ m" , n") = F(m', m" - m)F(n", n' - n)5 m+n i +n ^ n+m , +m „ , (17) 

where F(a,b) = (— l) a ^/(a + b)\/a\b\. In Fig. EI we illustrate the restrictions on 
arguments of the kernel (Tl6|) under which it can take non-zero values. 

Let us now consider how the nontrivial relation ( floT) looks like for the kernel ( fl6l) . 
Substituting (TTTT) for K in ( fl5l) . we obtain the equality 

F (k, m'" - m)F(m', m" - k)F{n"\ k + w!" - m - n'")F(n", w! + m" - n" - k) 

k 

= F ( m ', l + n'-n- m')F(l, n' - n)F{m", n" + n'" - m" - l)F(n"', n" - I) (18) 
i 

which is to hold true whenever m! + m" + m'" + n = n' + n" + n'" + m. Expressing n' 
through the other variables and using the explicit form of the function F yields 

/n"4-T)"'— m" n"\ 

{ n +nm ,n | (-!)'(/ + m > + m " + m '» - n " - n'" - m)\{n" + Tl'" - /)! 

^ + m" + m'" - n" - ri" - m)\(n" + ri" - m" - l)\(n" - l)\ 

l=\0,n"+n"'+m-m"-m"'] V ' K 
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(-l) n " +m "m"\(m' + m" + w!" - n" - n'" - m)\ 
im!" — m)\n"\ 

{m''m'+m"—n"\ , . , ., ,,, . ,. , ,, 

v 1 ^ (-1) (/c + m - m)\{m' + m" - k)\ 

X ^ k\(k - m - ri" - m"')\(m" - k)Um' + m" - n" - k)\ ' ' ' 

k=\0,m+n'"+m'"~\ K / v / \ / 

where \a,b] = max(a, b) and {a, b} = min(a, b). Although the obtained relation is 
rather complicated, it is valid for all non-negative integers m, m', m", m'", n, n", n'" . 
For instance, if we put m' + m" — n" = 0, then the summation over k in the right-hand 
side of Eq. ( Fl9l) reduces to a single term k = provided m + n'" — m'" < 0. As a result 
we derive a new property 

{n " ' V ^ {-l) l {M + l)\{n" + n"'-l)\ = , 

^ l\(M-m , + l)\(n' ,, + m'-l)\(n ,, -l)\ 1 ' ' 1 ' 

l=\0,m'-M] K M 7 v ' 

where M = m'" — m — n'" > and n" > m'. Note that the result of summation does 
not depend on m",n",n"'. In particular, if we choose n'" = 0, then 

y (-D'Cm + qi = ( _i r ' (2D 

l=lO,m'-M] K ' K ' 

Similarly, if we fix m' = in Eq. (I20I) . then 

{ v" } (-i)V + " w -0' = i r22 ^ 

^ i!(n w -Z)!(n w -Z)! ' 1 ' 

3.1. Intertwining kernel between star-product schemes 

Along with the considered star-product scheme with dequantizers ([6]) and quantizers 
©, there exist many other star-product schemes. In order to distinguish them let us 
use superscripts (91) for the just developed scheme of the normally ordered moments 
and write UW (m,n) and D^(m,n). 

As an example of the other quantization on the lattice (m, n) we may consider 
the star-product scheme with identical dequantizers and quantizers U^'(m',n') = 

(mf , n') = \m')(n'\. The symbols, that correspond to the two different schemes, 
are related by virtue of formulas 

oo 

ff\m,n)= K^{m,n-m',n')ff{m\n'), (23) 

m',n'=0 
oo 

ff{m', n') = n'; m, n)ff \m, n), (24) 

m,n=0 

where the intertwining kernels are expressed through dequantizers and quantizers as 
follows: 



K^ m (m,n;m',n') = Tr[U™(m 7 n)D®\m',n')] = ( ^^, <W,n+„', (25) 

K m ^(m',n';m,n) = Tr[U^(m' ,n')D^\m' ,n')} = ^ < WmW (26) 

(m — n')W m'\n'\ 
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It is worth mentioning that the dequantizers and quantizers of a star-product 
scheme should not necessarily depend on discrete variables. For example, in optical 
tomography, the dequantizer is a projection on the rotated quadrature, i.e. U^°\X, 0) = 
\X, 9) (X, 9\, where (gcos^ + psm9)\X, 6) = X\X,8), q and p are the position 
and momentum operator, respectively. The intertwining formulas that connect the 
optical tomogram and the normally (antinormally) ordered moments are derived in the 
papers [221I3Q1I35]. 



4. Evolution equations 

As it is outlined in Section (TJ the star-product formalism is quite useful for describing 
evolution of quantum states. In case of the star-product scheme ©-(111), one deals with 
the functions of discrete variables instead of operators. The more important fact is 
that the functions f p (m,n;t) are experimentally measurable [32H34] at different time 
moments t, giving an opportunity to observe dynamics of the system and motivating us 
to derive the evolution equations in terms of measurable quantities. 

To start with, the time dependent and stationary von Neumann equations for the 
density operator p take the following form within the star-product formalism: 

9 -l = -i[H,p) O dfp{m dt H] t] = -i {f H */,-/,* f H ) (m, n; t), (27) 

^ [Hp E + p E H} = Ep E & - (f H * f PE + f PE -k f H ) (m, n) = Ef PE (m, n), (28) 

where H and E are the Hamiltonian and the permitted energy level, respectively; the 
Planck constant h~ = 1. 



4-1. Moments' dynamics for harmonic oscillator 

Let us consider a free evolution of the electromagnetic field governed by a harmonic 
oscillator Hamiltonian H = a) a + |, where we have also used dimensionless units for 
the frequency (u> = 1). In this case we immediately encounter a problem of finding the 
symbol /i?(m, n) because it takes infinite values if m — n. However, this difficulty can 
be avoided if we focus on the star product of symbols. In fact, 

(fa * f P ) (m, n) = Tr[(a t ) m a n (a t a + \)p] = Tr[(tf) m (a n tf)ap] + \Ti[(a)) m a n p] 
= Tr[(a t ) m (a t a n + na n - l )ap\ + \Ti[{a)) m a n p\ 

= f p {m + 1, n + 1) + (n + \)f P {m, n). (29) 

Similarly, we obtain (f p ★ fn) (m, n) = f p (m + 1, n + 1) + (m + ^)f p (m,n). 
Substituting these results in Eqs. (|2"T )) - (|2"8~]) yields 

df p (m,n;t) 

— = i(m-n)f p (m,n;t), (30) 

Tfi -\- Tl ~~\~ 1 

f PE (m + l,n + l) + f PE (m,n) = Ef PE (m,n). (31) 
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Figure 3. Local dynamics of the moment {{<v) m a n ) on the "phase-space" lattice: 
(a) harmonic oscillator and (b) damped harmonic oscillator. Squares denote nodes for 
which the corresponding moments are involved in the evolution of ((a^) m a n ). 



It is not hard to check that symbols f^^^m^n; t) = N\e i{m - n)t 5 mtr J{N - m)\ of 
the Fock states \N) such that m = n < N do satisfy the derived equations (|30l)-(l3TT) if 
E = N+\. 

In case of the harmonic oscillator Hamiltonian, the dynamics (1301) of quantum states 
on the lattice (m, n) reduces to f p (m, n; t) = f p (m, n; o)e^ m_n ^, i.e. the moments simply 
gain phases in accordance with their position on the lattice. The corresponding vectors 
in Fig. [lb save their length and rotate with frequencies u mn = m — n, so in time interval 
t = 2tt all the vectors come back to the initial position. Thus, evolution of the moment 
({(v) m a n ) is extremely local on the "phase-space" lattice and does not depend on the 
values of moments in surrounding nodes (Fig. [3^). 

4-2. Moments' dynamics for damped harmonic oscillator 

Let us consider a damped evolution equation for the density operator described by the 
usual master equation (see, e.g., [53|[5^]) 



— = —i[d^d, p] + 7(1 + z/)(2dpa^ — a)ap — pa' a) + 7z/(2d^pd — aa) p — pdd^), (32) 

L/ v 

where 7 is the damping coefficient and v is the equilibrium mean number of photons in 
a given mode. 

Arguing as above, we find the star product of symbols in question 

(fa * f P * fat) {m, n) = f p (m + 1, n + 1), (33) 
(/at * f P * fa) (m, n) = f p (m + 1, n + 1) + (m + n + l)f p (m, n) 

+ mnfpirn — l,n — 1). (34) 

Now one can rewrite fl32l) in terms of the measurable moments ((d^) m d n ) = f p (m, n) 
as follows: 

df P (m,n;t) 



dt 



[i(m — n) — 7(777, + n)\ f p (m, n; t) + 2 r yvmnf p (rn — 1, n — 1; t). (35) 
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Figure 4. Nodes on the "phase-space" lattice for which the corresponding moments 
are involved in the evolution of a particle: (a) kinetic term p 2 /2, V(q) oc q 2 ; (b) 
V(q) oc q 3 ; (c) V(q) oc q 4 . The lower sum m + n for the marked nodes the greater 
weight of corresponding moments in the time evolution. 



In case of the damped harmonic oscillator, the dynamics of moments is again quite 
local on the lattice (m, n) and involves only two nodes (Fig. [3b). It is worth mentioning, 
that Eq. ( 1331) could also be derived by using a connection between the normally ordered 
moments and the Wigner function, and then substituting these relations in the Fokker- 
Planck equation for the Wigner function [35] . However, as we can see, the star-product 
approach is straightforward to derive evolution equations for measurable quantities like 
the normally ordered moments of the creation and annihilation photon operators. 

4-3. Moments' dynamics for a particle 

Although the background of our consideration is measuring radiation fields at 
microwaves, the developed star-product formalism can be also applied to the one- 
dimensional motion of particles governed by the Hamiltonian H = p 2 /2 + V(q) = 
— la — a' ) 2 /4 + V((a + tf)/y/2). The main point is that the time evolution of moments 
remains local on the lattice if the potential energy V(q) can be approximated by lower- 
order terms of the Taylor expansion (Fig. HJ). 

If V(q) = q , then the nodes (m',n') involved in the dynamics of the moment 
((a^) m a n ) form a truncated square with center at the point (m — 1, n — 1). The nodes 
(m', n') satisfy the relations 



m! — m + 11 + \n — n + 1 




. . . , / for even /, 
. . . , I for odd I, 



max(m' — m,n' — n) > 0. 

(36) 

It is instructive to compare the ^-function evolution and the evolution of moments. 
The derivatives d l V/dq l and powers q l form the Taylor expansion of V(q). These 
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derivatives show which partial derivatives d l W/dp l contribute to the evolution on the 
conventional phase space (q,p), whereas the powers q l show which moments ((a^) m a n ) 
contribute to the evolution on the "phase-space" lattice (m, n). The important practical 
difference between a quasidistribution on the conventional phase space (q,p) and our 
proposal is that the time development of a conventional quasidistribution is given in 
terms of the infinite-order partial differential equation which is not so easy to solve. A 
numerical solution of such an equation would require a construction of a two-dimensional 
rectangular grid in the (g,p)-plane, with the size and density of the grid being taken 
according to a desired accuracy (see, e.g., [IP]). The partial derivatives are then replaced 
by finite differences. The higher the order of the derivative, the more nodes of the grid 
are involved. Dealing with the ordered moments, one does not need to introduce any 
artificial grid because of the lattice itself. The time development is then given by the 
exact difference equations in contrast to the approximate finite-difference equations for 
quasidistributions on the (g,p)-plane. 

5. Conclusions 

To conclude, we present the main results of the paper. 

An analysis of the star-product scheme based on the normally ordered creation 
and annihilation photon operators has been motivated by the recent advances in 
measuring ordered moments for microwave quantum states [32TI34"] . In addition, the 
phase space of the constructed quantization scheme is a two-dimensional lattice (m, n) 
whose nodes are given by two nonnegative integers. Such a structure of the phase 
space is advantageous for describing the time development of some quantum systems 
because the exact evolution equations take the form of difference equations in contrast to 
the partial differential equations for conventional quasidistributions on the (g,p)-plane 
usually approximated by finite-difference equations on the (g^, p 3 -)-grid of rather artificial 
size and density. 

Moreover, it is quite reasonable to define a quantum state evolution in terms of 
the measurable quantities, so we have filled a gap of such equations in terms of the 
normally ordered moments ((a)) m a n ). The dynamics of moments is shown to be local 
on the lattice for (damped) radiation fields and particles moving in smooth potentials. 

Another substantial result is that the star-product kernel is found in the explicit 
form (Tl6|) . As any star-product kernel is to satisfy specific non-linear equalities, we have 
applied one of those equalities to the obtained kernel and thus derived new sum relations 
involving factorials ( fX9|) — ( |22]) . 
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